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Abstract: In this work we estimate the heavy quark production in the interaction of 
ultra high energy cosmic rays in the atmosphere, considering that the primary cosmic ray ^ ^ 

is a proton or a photon. At these energies the saturation momentum Q^ati^) stays above 
the hard scale = 4m^, implying charm production probing the saturation regime. In 
particular, we show that the ep HERA data presents a scaling on Tc = (Q^ + ^'^) / Ql^^{x) . o 
We derive our results considering the dipole picture and the Color Glass Condensate for- C !) 

malism, which one shows to be able to describe the heavy quark production in 7p and pp 
collisions. Nuclear effects are considered for the scattering of primaries with the air nuclei 
and we provide a parametrization for the charm and bottom differential cross sections, 
da/dxp, which can be used as an input for numerical implementations for lepton flux. Im- 
plications on the flux of prompt leptons at the Earth are analyzed and a large suppression 
is predicted. 
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1. Introduction 

Ultra high energy cosmic rays (UHECRs) remain a puzzle in physics. Although the exis- 
tence of UHECRs with energies above lO^o eV is now a well-established fact, the theoretical 
understanding of its origin and propagation is a subject of strong interest and intense dis- 
cussion Q. In general terms the current theories for the origin of UHECRs can be classified 
in two distinct scenarios. In the Bottom-Up scenario, charged particles are accelerated from 
lower energies to high energies in astrophysical environments. On the other hand, in the 
Top-Down scenario, the energetic particles arise from decay of massive particles originat- 
ing from physical processes in the early Universe. Another open question is the basic 
composition of the UHECRs. While in the Bottom-Up models the primaries are acceler- 
ated protons (or nuclei), the Top-Down models predict an increasing photon component 
at energies above 10^^ '' eV. Indeed, even in bottom-up models ultra high energy photons 
are expected from the GZK process during the cosmic ray propagation [||. Finally, the 
interaction of the UHECRs with the atmosphere nuclei probes the theory of the strong 
interactions in a new kinematical range characterized by a center of mass energy of ap- 
proximately 500 TeV, which is more than one order of magnitude larger than the future 
Large Hadron Collider at CERN. 

An important subject in cosmic ray physics is the flux of prompt leptons at the Earth 
which reflects the primary interactions at energies that can by far exceed the highest 
available energies. Its quantification is essential for the cosmic ray physics as well as for 
neutrinos physics. On one hand, the fiux of cosmic ray muons in the atmosphere, under- 
ground and underwater provides a way to testing the inputs of nuclear cascade models, that 
is, parameters of the primary cosmic ray fiux (energy spectrum, chemical composition, ...) 
and particle interactions at high energies On the other hand, the flux of atmospheric 
neutrinos and muons at very high energies provides the main background of searches for 
the muons neutrinos from extra-galactic neutrino sources in the neutrinos experiments 



(e.g. AMANDA, Antares, Nestor), limiting the sensitivity of any neutrino telescope to 
astrophysical signals Q. 

The flux of prompt leptons is directly associated to the charmed particle production 
and its decays, whose estimation is strongly dependent on the model used to calculate the 
charm production cross section and energy spectra 0, It is related to the need of 
extrapolating charm production data obtained at accelerators energies to the order of mag- 
nitudes higher energies of the relevant cosmic rays collisions. Another uncertainty present 
in the estimate of prompt lepton fluxes is associated to the fact that different authors do 
not use the same atmospheric particle showering routines, turning the comparison among 
their predictions even more difficult (See e.g. Ref. Q). 

Over the past few years, the heavy quark production at colliders has received consid- 
erable attention. In particular, the compatibility among perturbative QCD (pQCD) and 
experimental data was verified and today one can conclude that even for charm, whose 
low mass and larger non-perturbative hadronisation corrections might cast doubts on the 
predictive power of the theory, pQCD actually manages to deliver good results. This 
framework was used in the prompt lepton fiuxes calculations presented in Refs. ||5|, 0, 
However, since the lepton fiuxes are strongly dependent on the behavior of the gluon dis- 
tribution at small-x [10| and it is determined by the QCD dynamics, we have that any new 
dynamical effect will modify the estimates of the lepton fiuxes. Recent results at the DESY 
electron-proton collider HERA indicate the presence of the parton saturation effects (See 



e.g. Ref. [Q), which modify the linear DGLAP dynamics [11|. Furthermore, a current 
open question is related to the possibility of the breakdown of the collinear factorization 
at higher energies due to saturation effects which are expected to be present in this regime 
fl^ , 14 1 . As in the previous calculations ^, |7|, ^ the authors have assumed the validity 



of the collinear factorization and gluon distributions which are solution from the DGLAP 
evolution equation, we believe that a new study of the prompt lepton fiux is timely and 
necessary. 

In this paper we estimate the heavy quark production considering the present under- 
standing of the high energy regime of the theory of strong interactions (For recent reviews 
see Ref. [12 1). In this regime, perturbative Quantum Chromodynamics (pQCD) predicts 
that the small-x gluons in a hadron wavefunction should form a Color Glass Condensate 
(CGC), which is characterized by the limitation on the maximum phase-space parton den- 
sity that can be reached in the hadron wavefunction (parton saturation) , with the transition 
being specified by a typical scale, which is energy dependent and is called saturation scale 
Qsat- In order to estimate the saturation effects we calculate the heavy quark production 
using the color dipole approach, which gives a simple unified picture for this process in 
photon-hadron and hadron-hadron interactions. Beside charm production we also estimate 
the bottom production, since S-hadrons can also contribute to the r and u-r fiuxes. It is 
important to emphasize that differently of the charm quark production in the current accel- 
erators, where the saturation scale Qsat is smaller than the typical hard scale, ^iq = 2mQ, 
at the energies of interest in this paper, E > 10^*^ eV, we will probe for the first time 
the kinematical regime where Qsat > ^c- Therefore, at these energies one can expect a 
large modification of the charm quark total cross sections and, consequently, on the fiux of 



prompt leptons at the Earth. 

Here we restrict ourselves to estimate the heavy quark production considering the 
present understanding of the high energy regime of the theory of strong interactions and 
to predict the magnitude of the saturation effects in the main quantities which are used 
as input in the atmospheric particle shower routines. As the theoretical predictions of the 
prompt leptons depend strongly on the behavior of the total cross section for heavy quark 
production at high energies, we believe that this partial calculation allow us to obtain a 
reasonable estimate of the magnitude of the saturation effects on prompt lepton flux at 
the Earth. Of course, more detailed studies are necessary in future in order to get precise 
predictions. This paper is organized as follows. In the next section we present a brief 
review of the heavy quark production in the color dipole picture, demonstrating the direct 
relation between the total heavy quark cross sections and the dipole-target cross section 
adip- The QCD dynamics is discussed in Section ^ and the phenomenological model for 
adip used in our calculations is presented. In the Section ^ we present our results and our 
main conclusions are summarized in the Section |5[ 



2. Heavy quark production in the dipole picture. 

Let us start the analyzes considering the heavy quark production in photon-hadron inter- 
actions at high energies. It is usually described in the infinite momentum frame of the 
hadron in terms of the scattering of the photon off a sea quark, which is typically emitted 
by the small-x gluons in the proton. However, in order to disentangle the small-x dynamics 
of the hadron wavefunction, it is more adequate to consider the photon-hadron scattering 
in the dipole frame, in which most of the energy is carried by the hadron, while the photon 
has just enough energy to dissociate into a quark-antiquark pair before the scattering. In 
this representation the probing projectile fluctuates into a quark-antiquark pair (a dipole) 
with transverse separation r long after the interaction, which then scatters off the proton 



1 15]. In this approach the heavy quark photoproduction cross section reads as, 

aijp^QQX)= j dzcfr\^liz,r,Q^ = 0)\^adipix,r), (2.1) 
where '^j^iz, r, Q'^) is the transverse light-cone wavefunction of the photon which is given 

by 

1^^ {z, r, Q2 = 0)p = ^ el {[z' + (1 - zf] Kfier) + K^.ier)} (2.2) 

where the quantity = z{l — z) Q'^ + mg depends on the heavy quark mass, mq. The 
i^o,i are the McDonald function. The variable r defines the relative transverse separation 
of the pair (dipole), whereas z and z = (1 — z) is the longitudinal momentum fractions 
of the quark and antiquark, respectively. The dipole cross section, adip, parameterizes the 
cross section for the dipole-target (nucleon or nucleus) interaction. As usual, the Bjorken 
variable is denoted by x. 

Concerning heavy quark hadroproduction, it has been usually described considering 
the collinear factorization, with the behavior of the total cross section at high energies being 



determined by the gluon-gluon fusion mechanism and the gluon distribution at small- a;. 
However, at large energies this process can be described in terms of the color dipole cross 
section, similarly to photon-hadron interactions, using the color dipole formalism ||l^. In 
this approach the total heavy quark production cross section is given by 111 



/-ln(H^) 

a{pp ^ QQX) =2 ' dyxiG x a{GN ^ QQX) , (2.3) 

JO 

where y = ^ln(xi/x2) is the rapidity of the pair, ^ip ~ uiq is the factorization scale, 
x\G{xi, is the projectile gluon density at scale and the partonic cross section 
cj{GN QQX) is given by 

a{GN^QQX) = J dz d^r \^a^QQ{z,r)f aggG{z,r), (2.4) 

with "^Q^QQ being the pQCD calculated distribution amplitude , which describes the de- 
pendence of the \QQ) Fock component on transverse separation and fractional momentum. 
It is given by 

\^G^QQi^,r)\' = {mlKlimQT) + [z^ + z^] m^K?(mQr)} , 

where asipfi) is the strong coupling constant, which is probed at a renormalization scale 
fifi ~ niQ. Moreover, aqqc is the cross section for scattering a color neutral quark- 
antiquark-gluon system on the target and is directly related with the dipole cross section 
as follows 

9 1 

O-qqG = o l^dip{x2, Zv) + (Tdip{x2, Zv)] - -adip{x2,r). (2.5) 
o o 

The basic idea on this approach is that at high energies a gluon G from the projectile hadron 
can develop a fluctuation which contains a heavy quark pair [QQ). Interaction with the 
color field of the target then may release these heavy quarks. In Ref. |17] the equivalence 
between this approach and the gluon-gluon fusion mechanism of the conventional collinear 
factorization has been demonstrated. In particular, the dipole predictions are similar to 
those obtained using the next-to-leading order parton model calculation. Moreover, the 
predictions for heavy quark production agree well with the current experimental data. 

It is important to emphasize that the color dipole picture is valid for high energies, 
where the coherence length l^. ~ l/x2 is larger than the target radius. At smaller energies, 
we should include, for instance, the quark-antiquark contributions. However, at high ener- 
gies the dipoles with fixed transverse separations are the eigenstates of interaction in QCD, 
which becomes the eikonalization an exact procedure, allowing a direct generalization of the 
approach to calculate the heavy quark production in photon-nucleus and proton-nucleus 
collisions. Moreover, this approach is the natural framework to include the saturation ef- 
fects. Another aspects that should be emphasized are that (i) in this formalism a X-factor 
is not necessary, since is includes all higher- twist and next-to- leading corrections; (ii) the 
dipole cross section is universal, i.e., it is process independent. 



3. High Energy QCD Dynamics. 



The Color Glass Condensate is described by an infinite hierarchy of coupled evolution 
equations for the correlators of Wilson lines In the absence of correlations, the first 
equation in the Balitsky-JIMWLK hierarchy decouples and is then equivalent to the equa- 
tion derived independently by Kovchegov within the dipole formalism In the CGC 
formalism a^ip can be computed in the eikonal approximation, resulting 



adip{x, r) = 2 j cPhM{x, r, b) , (3.1) 



where M is the dipole-target forward scattering amplitude for a given impact parameter b 
which encodes all the information about the hadronic scattering, and thus about the non- 
linear and quantum effects in the hadron wave function. It is useful to assume that the 
impact parameter dependence of M can be factorized as AA(x, r, b) = Af{x, r)S{b), so that 
o'dip{x,r) = aQj\f{x,r), with (Tq being a free parameter related to the non-perturbative 
QCD physics. The Balitsky-JIMWLK equation describes the energy evolution of the 
dipole-target scattering amplitude J\f{x,r). Although a complete analytical solution is 
still lacking, its main properties are known: (a) for the interaction of a small dipole 
(r <^ 1/Qsat)i AA(r) r^, implying that this system is weakly interacting; (b) for a 
large dipole (r S> l/Qsat), the system is strongly absorbed and therefore M{r) « 1. The 
typical momentum scale, Qg^t ^ a;"^ (A ~ 0.3), is the so called saturation scale. This 
property is associated to the large density of saturated gluons in the hadron wave function. 
In our analysis we will consider the phenomenological saturation model proposed in Ref. 
||20|] (including fit with charm quark) which encodes the main properties of the saturation 
approaches, with the dipole cross section parameterized as follows 



AAo(^)'"^"^\ forf<2. 



a^^^{x,r) = ao 

1 - exp [-a In^ (6r)J , for r > 2 , 



where f = rQsat{x) and the expression for f > 2 (saturation region) has the correct 
functional form, as obtained from the CGC formalism ||l8[. For the color transparency 
region near saturation border (r < 2), the behavior is driven by the effective anomalous 
dimension 7efr (x, r) = 7sat + Jxi^l^i/x) ' "^here 7sat = 0.63 is the LO BFKL anomalous 
dimension at saturation limit and k = 9.9. Hereafter, we label this model by CGC. 

In the case of photon-nucleus and hadron-nucleus interactions we take the Glauber- 
Gribov formalism for nuclear shadowing, with an average mass number A = 14.5. This is 
done by replacing a dip for a nucleon by a nuclear dipole cross section. 



nucl o / j2t 



j d?b I - exp (^]^ATA{h)(Tdi^ 



(3.2) 



where Tji{b) is the nuclear profile function. This approach leads a reasonable description 



of the current experimental data for the nuclear structure function |21]. 

As a short comment, the present work is based extensively on the Color Glass Con- 
densate. However, there are complementary approaches in literature. For instance, in Ref. 



1 22] a loop-loop correlation model is able to compute hadronic cross sections in good agree- 
ment with cosmic ray data from Fly's eye and Akeno. Moreover, saturation in the impact 
parameter space can be explicitly investigated in that model. 
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Figure 1: Experimental data on inclusive charm production at DESY-HERA plotted versus the 
scaling variable t^. 



4. Results for Cosmic Rays 



Before presenting our results for the heavy quark production at very high energies let us 
consider a basic property of the saturation physics: the geometric scaling. When applied 
for the deep inelastic scattering it means that the total "f*p cross section at large energies 
is not a function of the two independent variables x and Q, but is rather a function of the 
single variable r = Q"^ /Ql^^{x) as demonstrated in Ref. [^]. In Ref. |24| we demonstrate 



that this property is also present in the charm experimental data at Q'^ > 0. In order to 
include the experimental data for the photoproduction of charm we have reanalyzed our 
results and verified that the DESY-HERA experimental data |2^ at > presents the 
geometry scaling behavior on the variable = (Q^ + Am'^)/Q'^^^. It is demonstrated in 
the Fig. where we also present the predictions of two scaling models |2^], which we 



generalize for charm production. In the symmetric saturation model |27], scaling on Tc has 
been computed before in Ref. ||2^ and here it is generalized for the case where saturation 
scale is larger than It reads as, 



crrot (re) 



1 — exp 



(l + log(Te) 



for Tc > 1 , 



i - :^ exp [-kTc (1 - log (tc)] , for Tc < 1 



The parameters are taken from [27|, with k = 0.8, and the overall normalization gives a = 
28 fih. Furthermore, we consider the model proposed in Ref. [^] (hereafter ASW model), 



which relate the ^p{A), pA and AA cohisions through the geometric scaUng property of 
the saturation physics. The parameters for this model has been obtained from a fit to the 



small-x ep DESY-HERA data and lepton-hadron data [28|. When generalized for charm 
production we obtain 



CT. 



7* A 



TTRj 



7i? + r 0, ^ + In 



'A,c, 



(4.1) 



where 7^ is the Euler constant, F (0, (3) the incomplete Gamma function and Ra is the 
nuclear radius. The nuclear scaling variable for charm production is ta.c = Tc [ ) • 
The parameters are a = 1.868, b = 0.746, A = 1.266, ttR^ = 1.55 fm^ and ao = 20.28 
//b. Our predictions for charm production in jp processes are presented in Fig. 0. The 
agreement with experimental data is very good, which allow us to propose the Eq. Q 
as a theoretical parameterization for the charm production in lepton-nuclei interactions at 
very high energies. 
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Figure 2: Profile function for charm production. 

We are ready now to investigate the charm production initiated by primary cosmic rays 
photons or protons interacting with atmosphere. Initially, lets consider the photon-hadron 
interaction and analyze the behavior of the profile function, defined by 

1 daQQ _ 47rr f ,2, , ,,.,7/. . ^2 



crtot dr 



^ j d^bdz\-i'l{z,r,Q^ = 0)\^ 1-exp (^-^ATA{b) a^ipix, r^)^ 



(4.2) 



which allow us to estimate the mean dipole size dominating the nuclear heavy quark pho- 
toproduction. In Figs. § and ^ are shown, respectively, the overlap function for the charm 
and bottom production as a function of dipole size. They are computed for the average 
mass number A = 14.5 and for different photon energies. In the charm case, the distri- 
butions are peaked at approximately r ~ 0.4 GeV~^, whereas for the bottom case this 



value is shifted to r ~ 0.1 GeV~^, which agree with the theoretical expectation that the 
gg pairs have a typical transverse size ~ l/2mj [15|. Therefore, the main contribution to 
the cross section comes from the small dipole sizes, i. e. from the perturbative regime. In 
contrast, for light quarks a broader r distribution is obtained, peaked for large values of the 
pair separation, implying that nonperturbative contributions cannot be disregarded in that 
case. We also show in the figures the corresponding nuclear saturation radius, estimated 
in a simple way as iisat — /Qsa.t{x), where x = -^^^ and W^j^ = 2m]^E^, for two 

different values of the photon energy E^. For the charm case we have that at low energies 
{Ery = 10'* GeV) the distribution is peaked at values smaller than the corresponding i?sat; 
which implies that the saturation effects can be disregarded in this regime. On the other 
hand, for high photon energies {E^ = 10^^ GeV) we have that the peak occurs at larger 
values than i?sat and the saturation effects are expected to contribute for the total cross 
section. In contrast, for the bottom case, we have that for all energy range considered the 
peaks occur at smaller values than -Rsati which implies that these effects do not contribute 
too much for the bottom contribution. 
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Figure 3: Profile function for bottom production. 



In Figs. ^ and ^ we show the charm and bottom photoproduction cross sections as a 
function of photon energy, respectively. The 7p interaction is obtained via Eq. ( |2.lD with 
the dipole-proton cross section given by c^jp*^- The experimental data from HERA |25| 
and fixed target collisions |26| are also included for sake of comparison. The theoretical 
curves roughly describe the data, which have large uncertainties and large errors. We call 
attention that the CGC model successfully describes the charm and bottom electroproduc- 
tion data (even at small Q^) and so we believe in its robustness. Therefore, accurate new 
experimental measurements in this region is of particular importance. 

As explained in the previous section, for the photon-air interaction we take the Glauber- 
Gribov formalism for nuclear shadowing, with an average mass number A = 14.5 and adip 



-8- 



replaced by cr^i^^K Initially, lets discuss the nuclear charm photoproduction. In Fig. ^ the 
accelerator data |25] and the predictions of our theoretical parameterization, Eq. ([4.1|), 
are also presented. The energy dependence is mild at > 10^ GeV , giving a oc -E^'^^. 
The parameterization is numerically equivalent to the complete CGC dipole calculation, 
which can be useful as input in Monte Carlo implementations. In order to estimate the 
magnitude of the saturation corrections to the high energy process, we also show the result 
using the color transparency (long dashed curves) calculation, which corresponds to the 
leading twist contribution adip oc 1/r^). The deviation are about one order of 

mag nitude at 10^^ Qgy, increasing for higher energies. Therefore, in agreement with 

the previous analyzes of the profile function, parton saturation plays an important role in 
charm production by cosmic rays interactions when the primary particle is a photon. 
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Figure 4: Charm production in 7p and jA interactions. Data from Refs. 

Concerning to nuclear bottom photoproduction, our results are presented in Fig. |5|. 
We have that the contribution of the saturation effects is smaller, being of approximately 
10 % at very high photon energies {E^ ~ 10^^ GeV). This result agree with our previous 
analyzes of the profile function. An important aspect is that the growth with the energy 
of the bottom cross section is steeper than the charm case due to distinct physics which is 
dominant in the production process (linear versus saturation physics). It implies that the 
contribution for the prompt lepton fluxes associated to the bottom becomes more important 
for larger energies. 

Lets now consider that the primary cosmic ray is a proton and estimate the heavy quark 
production cross section in this case. In Fig ^ we show the theoretical predictions for charm 
quark production in proton-proton and proton-air interactions as a function of energy of 
primary cosmic rays. The solid curves correspond to the calculation using Eq. ( |2.3[ ), with 
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Figure 5: Bottom production in 7p and interactions. Data from Refs. |25: E6|. 



eluon distribution given by the NLO "physical gluon" MRST2004 parameterization pi. 
The accelerator data |Q and the color transparency (long dashed curves) calculation are 
also presented. We have that the current data at high energies are well described using the 
color dipole picture and that for the energy range of the data the color transparency and 
CGC prediction are almost identical. The energy dependence is suddenly steeper than the 
photon case, with a cx -Ep '^^ at Ep > 10^ GeV. The deviation between CGC and leading 
twist approaches remains large, being of order of six at Ep PS 10^2 and mcreasmg 

with energy. It confirms the importance of parton saturation corrections in charm quark 
production via cosmic rays. 



In Ref. 1 10] the authors have computed the total cross section for charm production 
using the collinear factorization, next-to-leading order corrections and different parton 
distribution parameterizations. In particular, the dependence of the predictions in the 
behavior of the gluon distribution at small-x was analyzed in detail. In comparison with 
those results we have that our predictions considering the color transparency limit are 
similar with the MRST one for A = 0.2 (See Fig. 1 in Ref. |jlO[). On the other hand, our 
predictions considering the saturation effects are very similar to those obtained previously 
by Thunman, Ingelman and Gondolo (TIG) ||5[, which have used an option of PYTHIA 
by which the gluon distribution is extrapolated for x < 10~^ with A = 0.08. As discussed 
before, the prompt lepton fluxes are strongly dependent on the charm cross section. Con- 
sequently, we can expect that if our results are used as input in the atmospheric particle 
shower routines, the predicted prompt lepton flux should be similar to the TIG predictions, 
which have been considered a lower bound for the flux of leptons obtained using pQCD. 
In other words, we have that the saturation physics on charm quark production implies a 
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Figure 6: Charm production in pp and pA interactions. Data from Ref. BC 



suppression of the prompt lepton fluxes. 
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Figure 7: Comparison of the CGC and color transparency predictions for the Xi?-dependence for 
charm production in pA interactions. 



In what follows we analyze the behavior of the XF-distribution for charm production, 
xpda^/dxF, which is necessary to calculate the charm production spectrum- weighted mo- 
ment, also called production "Z-moment". The Xi?-distribution can be directly obtained 



from the rapidity distribution in color dipole approach, which is given by |16] 



daipp -> QQX) ^ , ^ ^ 

dy 

Using that xp = xi — X2 and y = ^ln(xi/a;2) one obtain xp = — ^^sinhy, where Mqq 
is the invariant mass of the heavy quarks pair and the momentum fractions are given by 

xi^2 = (^Z + ^^^'^'^ ± xf)/2. Consequently, the Xi?-distribution is 

dajpp ^ QQX) ^ 1 dajpp ^ QQX) 

dxF . /^9A/r_ = /./^^2 . ,2 dy 



2Mqq/^sY + {xf) 

The extension to scattering on nuclei is straightforward by replacing a{GN — > QQX) for 
a{GA — > QQX). It is important to emphasize that the heavy quark production actu- 
ally contains two different phases, which are not unambiguously separable: the first is the 
description of the heavy quark in the hard collision which we have addressed above; the 
second is the non-perturbative fragmentation into a heavy hadron, which cannot be cal- 
culated and are in general extracted from the e+e~ data. In our study we do not include 
the hadronization of the heavy quarks, since it is model dependent. Of course, in a full 
calculation it should be considered. 

In Fig. 1^ we present this distribution for two values of the proton energy considering 
the saturation effects (CGC) and the color transparency limit. As expected, the difference 
between the predictions increases with the energy, being of approximately one order of 
magnitude for Ep = 10^^ GeV. In Fig. ^ we present the CGC predictions for the xp- 
distribution and different values of energy. The shape of the distribution becomes a little 
steeper as the energy increases. In order to allow future studies of the prompt lepton 
fluxes considering saturation physics we have parameterized the xjT'-distribution for charm 
production considering the following general form 

XF d^(P + Air^'^ + ^) =Ax<^Al- x^^r , (4.5) 

dXF 



with the quantities 



A = ao + 



ai In 



p 



108 GeV 



a = bo — bi In [ — ) , n = uq + rii In ' ^ 



104 GeV / ' ' ' V 104 GeV 



The parameters for two distinct energy ranges are presented in the Table ||. 

Let us consider now the bottom production cross section and its XF-dependence con- 
sidering the color dipole picture and saturation physics. From the previous discussion we 
have that these effects are expected to be minimal and that the cross section is dominated 
by the color transparency regime. However, we present our predictions obtained using color 
dipole picture and an eikonalized dipole-nuclei cross section, in order to establish the results 
which comes from this framework. As emphasized in Ref. the high energy production 
of u-r from beauty decays is not negligible, what motivates our study. In Fig. ^ we present 



Energy Range 


ao 


ai 


bo 


bi 


no 


m 


10^ <Ep< 10^ GeV 


6804 


826 


0.05 


0.016 


0.075 


-0.107 


10^ <Ep< IQii GeV 


5422 /ib 


403 


0.025 


0.023 


6.7 


-0.102 



Table 1: Parameters of the fit to the differential cross section xp da/dxp for charm production in 
proton-air collisions (see text). 
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Figure 8: Xi^-dependence for charm production in pA interactions predicted by saturation physics. 

our predictions for the energy dependence of the total bottom production cross section con- 
sidering the CGC and color transparency limit. As expected, we have that in the energy 
range considered in this work both predictions are almost identical for proton-proton and 
proton- nucleus interactions. In Fig. ^ the Xi?-distribution is shown for different proton 
energies. The behavior is similar to that obtained for charm production considering the 
color transparency limit, with the distribution becoming steeper for larger energies. As 
for the charm case we have parameterized this distribution considering the general form 
given in Eq. (|4.5| ) . The resulting parameters are presented in the Table |^ for two distinct 
energy ranges. It important to emphasize that similarly to the photoproduction case the 
bottom cross section increases with the energy faster than the charm one, due to distinct 
physics which is dominant in the production process. In order to illustrate this feature, we 
have estimated the ratio of the bottom and charm cross sections for proton-air collisions 
as a function of proton energy for Ep > GeV and verified that it is proportional to 
Rb/c Of Ep'^- Basically, we have that the bottom production cross section is increasingly 
important and reaches 10-13 % of the charm cross section for ultra-high energy cosmic 
rays. Consequently, the contribution for the prompt lepton fluxes associated to the bottom 
becomes important for larger energies. 




Figure 10: a;F-dependeiice for bottom production in pA interactions. 



5. Summary. 

The determination of the prompt lepton flux is fundamental in order to establish, for 
instance, the background for ultra high energy neutrinos from cosmological sources. At 
high energies the lepton fluxes are quite sensitive to the heavy quark cross section, which 
implies that the choice of an appropriate theoretical framework to estimate the cross section 
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Energy Range 
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Ol 


bo 


bi 


no 
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10^ <Ep< 10^ GeV 


180.23 /ib 


19.79 


-0.048 


0.035 


7.54 


-0.091 


10^ <Ep< 10" GeV 


236.13 /ib 


36.20 


-0.061 


0.023 


9.42 


-0.18 



Table 2: Parameters of the fit to the differential cross section xp da/dxp for bottom production in 
proton-air collisions (see text). 



in this energy range is fundamental. Previous calculations have considered perturbative 
QCD at next-to-leading order, assuming the validity of the colhnear factorization and of 
the DGLAP dynamics in the kinematical range of very high energies. However, current 
accelerator data already have indicated the presence of new dynamical effects, associated 
to saturation physics. As the contribution of these effects increases with the energy, we can 
expect that it cannot be disregarded in the description of the interaction of ultra high energy 
cosmic rays with the atmosphere. Here we have estimated the heavy quark production in 
the interaction of cosmic rays in the atmosphere taking the primary cosmic ray as a proton 
or a photon. At ultra high energies and charm production the saturation scale stays above 
the semihard scale ^'^ and the process contains sizable contribution from the saturation 
regime. In particular, geometric scaling for charm production on the scaling variable Tc 
is demonstrated using small-x DESY-HERA data. Within the color dipole approach and 
CGC formalism for dipole-target interaction, the parton saturation corrections are huge, 
suppressing the cross section by one order of magnitude for ultra-high energy primaries at 
~ 10^ GeV. We also predict a factor ten of suppression of the prompt lepton fluxes for 
a pure proton component for the primary cosmic ray, increasing if the primary component 
changes as appropriate for a Top-Down model. This suppression is also present in the 
x^-distribution, which one the main inputs to the calculate the prompt lepton fluxes. The 
resulting predictions for the total cross section at very high energies are similar to those 
obtained previously by Thunman, Ingelman and Gondolo, which have been considered a 
lower bound. As there is a strict relation between the charm production and the prompt 
lepton fluxes, we believe that the resulting lepton fluxes obtained using our predictions for 
charm production as input of the atmospheric particle showers routines should be similar. 
In other words, we expect a suppression of the prompt lepton fluxes associated to the 
saturation physics when compared with those resulting from NLO calculations using the 
collinear factorization. Of course, a more detailed analyzes is necessary in order to quantify 
precisely this suppression. It is important to emphasize that these predictions should be 
considered a lower bound in the suppression, since in the CGC formalism the breaking of 
the factorization for heavy quark production is predicted for Qsat ^ /^c |13l) which implies 
a larger suppression. 

We also estimate the bottom production considering the dipole picture and saturation 
physics. In this case, we have that the saturation effects can be disregarded and the color 
transparency limit determines the behavior of the total cross section in the energy range 
of interest in this paper. However, as the i?-hadron decays should contribute significantly 
for the fiux of high energy Uj- neutrinos, the quantification of this cross section using the 



color dipole picture, which is expressed in terms of the eigenstates of interaction in QCD, 
is useful. 

We present a parameterization for the x^-distribution for charm and bottom produc- 
tion resulting of our calculations. It can be useful for future calculations of the prompt 
lepton fluxes. The main conclusion of our phenomenological analyzes is that the saturation 
effects implies a suppression of prompt lepton fluxes. Of course, in a full calculation we 
should include the fragmentation of the heavy quark pairs into hadrons and their subse- 
quent semileptonic decays, as well as, the propagation of the high energy particles through 
the atmosphere. Wc postponed these analyzes for a future publication. 
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